Exact numerical simulations of a one-dimensional, trapped Bose gas 
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We analyze the ground-state and low-temperature properties of a one-dimensional Bose gas in a 
harmonic trapping potential using the numerical density matrix renormalization group. Calculations 
cover the whole range from the Bogoliubov limit of weak interactions to the Tonks-Girardeau limit. 
Local quantities such as density and local three-body correlations are calculated and shown to agree 
very well with analytic predictions within a local density approximation. The transition between 
temperature dominated to quantum dominated correlation is determined and it is shown that despite 
the presence of the harmonic trapping potential first-order correlations display over a large range 
the algebraic decay of a harmonic fluid with a Luttinger parameter determined by the density at 
the trap center. 
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Stimulated by the recent experimental progress in gen- 
erating ultracold trapped quantum gases in one dimen- 
sion 0, S 0, S S there is a growing interest in cor- 
relation properties of these systems. The physics of one- 
dimensional quantum gases is distinct from that in higher 
dimensions as it is dominated by quantum fluctuations. 
In a homogeneous system of Bosons there is no long-range 
order even at T — 0; correlations decay as a power-law 
due to zero-point phase fluctuations. At any finite T 
there is an asymptotic exponential decay. The most pe- 
culiar property of interacting Bosons in ID is the tran- 
sition to the fermion-like Tonks-Girardeau gas 0, § for 
small densities or large interactions. The transition is 
characterized by a single effective interaction parameter, 
the Tonks parameter 7, where small values correspond to 
the weak-interaction or Bogoliubov limit and large val- 
ues to the Tonks-Girardeau limit. The homogeneous gas 
is exact ly s olvable by Bethe ansatz for T = |9j and 
finite T |lOj . Correlation properties can however not eas- 
ily be extracted from the Lieb-Liniger solution and 
require in general numerical techniques such as Monte- 
Carlo simulations |l2| . Approximate analytic expressions 
can be obtained only for small distances |13| or within the 
harmonic-fluid approach 0, . 

In the presence of a trap potential V(x) integrability 
is lost. In order to nevertheless calculate local proper- 
ties Bethe-ansatz solutions for the homogeneous gas can 
be employed together with a local density approximation 
(LD A) and the Hellman-Feynman theorem 16]. Recently 
we have used stochastic simulation techniques to calcu- 
late the density profile and first-order correlations of a ID 
Bose gas in a harmonic trap 0] • The stochastic simula- 
tions were however limited to temperatures larger than 
the trap energy ^Tk Huj and thus did not allow to go 
deeper into the quantum regime. In the present paper we 
develop an alternative numerical approach based on the 
density-matrix renormalization group (DMRG) 0, 0] 
which leads to results with much higher precision for tem- 
peratures form zero to fko. 



We consider a one-dimensional Bose gas with delta in- 
teraction in a (harmonic) trapping potential V(x) 
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where we have used oscillator units, i.e. h = m = 1. g 
is the ID interaction strength proportional to the s-wave 
scattering length in one dimension oid H3- 

In the absence of an external trapping potential the 
Hamiltonian Q is integrable in the thermodynamic 
limit, i.e. it has an infinite number of constants of mo- 
tion. The ground-state solution for T = which can be 
obtained by Bethe ansatz |i| shows that the ID Bose-gas 
is fully characterized by one parameter 7 = g/p, the so- 
called Tonks parameter. Here p the density of the gas. 
The Bethe ansatz leads to the so called Lieb equation 
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where A is an implicit function of 7: A = 7 J _. dk a(k). 
All local properties of the gas can be expressed in terms 
of the (even) moments of cr(k) 
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dkk m a(k), m = 2,4, . 



E.g. the equation of state reads p — p(p, g) — g 2 f(j) 



/(7) 



3e 2 ( 7 ) _ 7^(7) 
2 7 2 



(2) 



Integrability is no longer given when a (harmonic) 
trapping potential V(x) is taken into account. An of- 
ten used approximation to nevertheless describe the local 
properties in the inhomogeneous case is the local den- 
sity approximation (LDA). The LDA assumes that the 
homogeneous solution holds with the chemical poten- 
tial p replaced by an effective, local chemical potential 
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fi e s(x) = jU — V(x). As long as the characteristic length 
of changes is small compared to the healing length the 
LDA is believed to work well. Within this approximation 
one finds e.g. for the density of the gas: 
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where / 1 is the inverse function of eq.(|2J). 

In order to develop an in principle exact numerical 
algorithm we here want to employ powerful real-space 
renormalization methods such as the DMRG HEEDl. To 
this end it is necessary to map the continuous to a lattice 
model. As shown in |l7j,|22j this can be done in a consis- 
tent way by introducing an equidistant grid Xj — jAx, 
j G Z, which amounts to replacing the field operator 
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where a, is a bosonic annihilation 



operator. Integrals are replaced by their corresponding 
sums and the second derivative in the kinetic energy term 
can be safely approximated by the difference quotient 
£ji>(xj) « (it(x j+1 ) + ttfo-x) - 24>(x 3 ))/Ax 2 . This 
leads to the Bosc-Hubbard Hamiltonian 



-J(oJfij_i + a\a i+ i) + A fit a j + yat 2 <r 
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where J = l/2Ax 2 , A = l/Ax 2 + V(xj) - fx and 
U = g/ Ax. Expressing the scaled hopping in terms of the 
Tonks parameter at the trap center, J/U = (7p(0)Ax) _1 
and taking into account that aiD -C Ax <C ^(O)^ 1 , the 
ID gas corresponds to a compressible phase of the BHM 
with negative effective chemical potential approaching 
the line (J-bh/U = —2 J/U. In the limit of vanishing aver- 
age particle number per site, pAx — > 0, Hamiltonian (JTJ 
and become equivalent. 

The numerical DMRG calculations of the density pro- 
file, shown in Fig. (I), for Tonks parameters 7 ranging 
from 0.4 to about 70 show excellent agreement with the 
Lieb-Liniger result with LDA J3J apart from a very small 
region at the edges and the barely visible Friedel-type 
oscillations, which result from the finite number of parti- 
cles. One recognizes the typical change of the density pro- 
file from an inverted parabola in the Bogoliubov regime 
7 -C 1 to the square root of a parabola in the Tonks- 
Girardeau limit 7>1 [2(J • 

An important consequence of the Fermion-like behav- 
ior of Bosons in the Tonks limit 7 ^> 1 is a dramatic 
reduction of the loss rate due to inelastic three-body col- 
lisions 0. The rate is proportional to the local three 
particle correlation gs{x) = {^ 3 {x)^ 3 {x)) / p{x) 3 , and 
determines the stability of the Bose gas. Making use of 
the Hellman-Feynman theorem and the constants of mo- 




FIG. 1: Density of the ID bosonic gas in a trap at T = 0. 
The solid lines are the DMRG results and the dashed lines are 
the Lieb-Liniger prediction in local density approximation. 
Increasing values of 7 correspond to decreasing densities at 
the trap center. Insert shows details at the edge of the density 
distribution. 



tion of the homogeneous gas Cheianov [25j has found 
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Fig. [3 shows a comparison between the numerical data 
for 33(0) at the trap center with eq.JHJ) and the asymp- 
totic expression in the Tonks- Girardeau limit with 7 
taken at the trap center 7(0) = g/p(0). One recognizes 
again excellent agreement except for a small deviation for 
very large 7, where the numerics is however susceptible 
to errors due to the smallness of 53 . 
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FIG. 2: Local third-order correlations as function of Tonks 
parameter at the trap center (red crosses) compared to pre- 
diction from Lieb-Liniger theory with local density approxi- 
mation (solid line) and Tonks- Girardeaux limit (dashed line). 

In contrast to local quantities, such as the mo- 
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merits of the number density, information about spa- 
tial correlations of the homogeneous ID Bose gas such 
as gi(x 1 ,x 2 ) = {^(x 1 )^{x 2 ))/y / p(x 1 )p(x 2 ) cannot 
straight-forwardly be obtained from the Lieb-Liniger 
and Yang- Yang theories. Making use of the Hcllmann- 
Feynman theorem and the asymptotic properties of the 
Lieb-Liniger wavefunction for large momenta, Olshanii 
and Dunjko derived the lowest-oder terms of the Taylor 
expansion of 51 around x\ — x 2 — 13] 



9i(x u x 2 ) = 1 - 2(^(7) -lt' 2 {l))p 2 x 2 
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with x = x\ — x 2 . In the presence of a trapping po- 
tential the Tonks parameter becomes space dependent 
7 — > "f(x). Short-range correlations are however expected 
not to depend on the confining potential. FigOl shows a 
comparison between g\ obtained from eq. 10 and numer- 
ical results for different Tonks parameters at the trap 
center. Taking into account that a high resolution of the 
short-distance behavior is numerically very difficult the 
agreement is rather good. 




0.8 

0.6 
0.5 



•--•N=109, y(0)=0.41 
mN= 50, 7(0)=1.35 
«--»N= 95, 7(0)=3.95 
Ar^N=101,7(0)=10.2 
< < N= 67, 7(0)=20.3 
>--»-N= 83, 7(0)=66.9 







0.05 



x/L n , 



0.1 



0.15 



FIG. 3: First order correlations (dashed lines) compared to 
analytic short-distance expansion (solid lines) for a homoge- 
neous gas with 7 taken at the trap center. Values of 7 increase 
from top to bottom curve. 

The long-range or low-momentum behavior of the cor- 
relations can be obtained from a quantum hydrodynamic 
approach [lil llol| in which long-wave properties of the 
ID fluid are described in terms of two conjugate vari- 
ables, the local density fluctuations 8p and the phase 

and d 



= v/ p(x)e 8< H X ). The equations of motion for 5p 
follow from the effective Hamiltonian [22], |23J 
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Here vn = (ir)~ 1 dp/dp and vj — irp. 

In the homogeneous case one finds that the leading- 
order term in the asymptotic of the first order correlation 



at temperature T is given by [22| 

K/L T 
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where K = \/vj/vn is the so-called Luttinger parameter 
and Lt is the thermal correlation length Lt — wp/KT, 
where we have set ks = 1- One recognizes that for T = 
correlations decay asymptotically as a power-law with 
exponent 1/2 A", while for finite T there is an intermediate 
power-law behavior turning into an exponential decay for 
|*Ei — x 2 \ > Lt- For T — the exponent 1/2K is given 
by 
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FIG. 4: First order correlations in the temperature regime be- 
tween exponential and algebraic decay, top: semi-logarithmic 
plot, bottom: double-logarithmic. Solid curves are DMRG 
calculations in the trap, dashed lines are harmonic fluid pre- 
dictions for a homogeneous gas with 7 taken at the trap cen- 
ter. Transition from thermal (exponential decay) to quantum 
dominated correlations (algebraic decay) at T -C u is appar- 
ent. The parameters are: 7 = 3.95, N = 12. Increasing 
temperature values correspond to lower curves for large x. 
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In Fig0] we have plotted the first-order coherence 
gx(x,—x) for symmetric positions with respect to the 
trap center for 7 = 3.95 and different temperatures. For 
comparison the harmonic-fluid results for the homoge- 
neous case, eq.©, are also shown with K and p taken at 
the trap center and for T = 0. (The change of K and p 
with T has little effect). One recognizes two things: First 
of all the transition from an exponential to a power-law 
decay happens around T = 0.1a;. Secondly the correla- 
tions are rather well described by the homogeneous so- 
lution (JSJ. A similar observation can be made at T = 0. 
Fig. [SJshows the DMRG results for g%(x, — x) for differ- 
ent interaction strength. The straight lines show the har- 
monic fluid predictions for the homogeneous case. Again 
a rather good agreement is found for x < 3L osc , which on 
first glance is rather surprising since the density of the 
gas is space dependent. The agreement is less surprising 
if one notes that the local Tonks parameter j(x) and thus 
the local Luttinger parameter K (x) are almost constant 
within this distance range. Furthermore replacing p in 
the denominator of eq.JHJ) by p(x) and expanding in a 
power series one finds that even in the Tonks limit where 
K — > 1 the corrections are small for positions sufficiently 
far away from the edges of the density distribution. 
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FIG. 5: Logarithmic plot of first-order correlations for T = 
and various interaction strengths (dots). The dashed lines 
show power-law prediction from the harmonic fluid approach 
with a Luttinger parameter determined by the density at the 
trap center. 7 increases from top to bottom curves. 

In summary we have developed a numerical scheme 
based on the density-matrix renormalization group that 
allows to calculate local properties as well as correlations 
of a ID Bose gas in a trapping potential for tempera- 
tures up to the oscillator frequency. For local quantities 
such as the density or the local three-body correlation 
we found excellent agreement with the predictions from 
the Lieb-Liniger and Yang- Yang theories with local den- 
sity approximation. Deviations from LDA are found only 
in the immediate vicinity of the edges of the gas or for 
smaller particle numbers where finite size effects come 



into play. We have shown that first-order correlations for 
positions well within the gas are well described by the 
homogeneous theory with parameters taken at the trap 
center. In particular the transition from a thermal dom- 
inated regime of exponential decay to a power law decay 
of correlations was shown, with exponents as predicted 
by the harmonic fluid approach in the homogeneous case 
for parameters taken at the trap center. 
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